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SUMMARY 

A disk-loaded circular wave guide structure and test fixture were fabricated. The dispersion 
characteristics were found by theoretical analysis, experimental testing, and computer simulation using 
the codes ARGUS and SOS. Interaction impedances were computed based on the corresponding 
dispersion characteristics. Finally, an equivalent circuit model for one period of the structure was 
chosen using equivalent circuit models for cylindrical wave guides of different radii. Optimum values 
for the discrete capacitors and inductors describing discontinuities between cylindrical wave guides were 
found using the computer code TOUCHSTONE. 


INTRODUCTION 

The disk-loaded circular wave guide structure is a forward-wave type structure and has 
applications in linear accelerators, high power TWTs, and millimeter-wave devices [1]. The theoretical 
analysis of such a structure can be accomplished rigorously by using a generalized mode matching 
technique. Scattering parameters for circular wave guide junctions have been calculated using a 
computer program based on this type of technique [2]. 

Another effective approach based on a TM 01 field analysis [3,4] has been done to determine 
useful parameters such as power flow, group velocity, and interaction impedance for this type of 
structure. One of the striking characteristics of these structures is the high gain per unit length, nearly 
twice that of a helix structure with the same beam hole size and the same beam velocity [3]. The disk- 
loaded circular wave guide structure, however, is limited in bandwidth when compared to the helix 
structure and would be most suited to high-gain, low power applications such as in a miniature 
traveling- wave device [3]. 


THEORETICAL EVALUATION 

The disk hole radius, the cylindrical wave guide radius, the slot length, and the disk thickness 
are represented by a, b, s and t, respectively (Figure 1). 
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Using cylindrical coordinates, the fields for k < 0 in the region p < a are represented by 


E z = E 0 3 0 ( 7 p)e ~ = E 0 l 0 (rp)e ~ (1) 

E p = jE Q {(3h)l x { 7 p)e ~ = jE Q {0/T)l l { T p)e ~ & (2) 

E* = 3 {EjTi)[kh)l x ( 7 p)e - = j{Ejri){k/r) l x {r p)t ~ (3) 

They are assumed to be finite, axially symmetric, and propagating in the 2-direction. The 
parameters 7, /?, and k are the radial, axial, and free space propagation constants with 
j 2 — k 2 - p 2 = — r 2 , and »; = (/i 0 /e 0 ) 1/1,2 = 377 fi. J Q and Jj represent zero and first order Bessel 
functions while I 0 and Ij represent zero and first order modified Bessel functions. 


In the region a < p <b we assume that the electric field is TEM p with the guide wavelength 
much greater than the period L such that E x = 0 at p = b. 

E x oc [J 0 (MN 0 (W) - N 0 (tp) J 0 (H)]e ' (4) 

E+ = 0, (5) 

cx Uh)[h( k P)M kb ) - Nj(*/»)J 0 (W)]« - i/3x , (6) 

Again, J 0 and J x represent zero and first order Bessel functions while N 0 and N x represent zero and first 
order Neumann functions, respectively. H ^ is found from V X E = — jup 0 H. 

Matching radial admittances {H^/E x ) across the slot length in both regions results in the 
dispersion relation for the disk-loaded structure as given by Chu and Hansen [4] 


! Jl(7«) lM™)_ 1 l J 1 (ta)N 0 (tt)-N 1 (ta)J 0 (tt) 

7a J 0 ( 7 a) ™ I 0 (ra) ~ 1 -t,ka 3 0 (ka)N 0 (kb) - N 0 (*a)J 0 (i6)' 1 > 

The term 1/(1 — »y) is an empirical correction to the zero disk thickness result. It involves the spatial 
factor t) = t/(t + s) which takes into account the fact that the electric field E x is zero at the disk edge. 
The electric field across the slot at p = a is assumed to be constant in the 2-direction. 

Taking the derivative from the above dispersion relation with respect to k we find 


where 


and 


d0 _k 7 da(k)/dk k . r da{k)/dk 

d k f} fi d<f>( 7 )/d 7 p P d<£(r)/dr 

(8) 

^ (7 ) "^J 0 (7 I 0 (ra) 

( 9 ) 

1 1 J a (fca)N 0 (W)-N 1 (fca)J 0 (fc6) 

{ ) -l-qka l 0 (ka)N 0 (kb) - N 0 (fca)J 0 (fc6)' 

(10) 


The group velocity can be obtained using dp/dk and is given by 




(/Vo) 


1/2 


1 

dp/dk' 


( 11 ) 
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The low frequency cutoff of the disk-loaded wave guide is simply that of the cylindrical wave 
guide without disks. It is determined from the condition J 0 (jfc6) = 0 which yields kb « 2.405. This 
represents a wavelength Aj c « 2x6/2.405. The high frequency cutoff of the structure depends on the 
disk hole size provided the period L<(b — a)/2. This cutoff frequency is determined from the 
condition J 1 (ta)N 0 (ifc) - N 1 (ia)J 0 (i6) = 0. The wavelength which corresponds to the high frequency 
cutoff is approximately given by A « 4(6 — a). If L> (b — a)/2 then the high frequency cutoff is 
determined solely by the period of the disks and A = 2 L. 

The appropriate expressions for the electric and magnetic field intensities, equations (2) and 
(3), are used to evaluate the time-average power flow through the disk-loaded structure. 

2ir a 

P = ±Xe{J(ExH*)-ds} = ^ | EpH^ p dp d <f> 

S 0 0 

p = * nj~ ^3T j h\rp)p dp = 7C I^ra) r • p d(rp) 

0 0 


P =~W~\^Z « 2 [ I l 2 ( ra ) - I o( ra ) I 2( ra )]- 


( 12 ) 


Using the expressions for E z and P, equations (1) and (12), respectively, the interaction 
impedance is found to be 

z{a) - \ E ^P)\ = n(r/l3)(r/k) 1 0 2 (rp) 

2/? 2 P *p 2 a 2 [I 1 2 (ra) - I 0 (ra)I 2 (ra)]" 


(13) 


The average interaction impedance is found to be 

/ \E x \p)\ dS 


Z(p)=- 


2 p 2 PS 


T)(T/f3)(T/k) l 0 2 (rp)-l 2 (rp) 

xf3 2 a 2 (l! 2 (ra) - I 0 (ra)I 2 (ra)]' 


(14) 


The area of integration is over the cross-sectional area S of the electron beam [5]. For the calculation 
of average interaction impedance, the cross-sectional area S = xr 2 was that of a perturbing dielectric 
rod of radius r = 0.0155 inch used in the experimental evaluation of the structure. 

Space harmonics for the axial field are ignored in this analysis but can be found for the z- 
component of the electric field. 

00 a 

. - 


E zM = A n l o( r n P) e 3 


(15) 


with t 2 = 0 n 2 — k 2 and f) n = f3 Q + 2 vn/L. Assuming a uniform electric field E a at p = a across a slot 
of width /, the coefficients A n can be found from 


E z (a,z) = A M r n a ) e 


-iK* 


= E a - 


(16) 


n = ~ oo 


Actually the electric field at p = a across the slot varies approximately as cosh(mz) where m is 
experimentally determined [6]. 
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or 


• • 1^3 j 

Multiplying by e m 


and integrating over the slot we have 



dz 


A m = E a 


l/L sin(/? m //2) MJ/ 2 ) 
Io( r ra a ) Pj! 2 


( 17 ) 


Equation (7) was evaluated to obtain the theoretical dispersion characteristics for structures with 
0.003, 0.005, 0.008, 0.010, and 0.020 inch thick disks and slot lengths of 0.16 inch (Figure 2). 



Figure 2 - Theoretical dispersion characteristics of the disk-loaded structures. 


Equation (14) was then evaluated to find the theoretical interaction impedances (Figure 3). 



0L/7T 

Figure 3 - Theoretical interaction impedances of the disk-loaded structures. 
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EXPERIMENTAL EVALUATION 


In many instances the theoretical analysis of a circuit may be too complicated to solve for the 
dispersion relation of a structure. The dispersion characteristics can be experimentally measured using 
a standard technique [1,7] in which resonance frequencies of a section of the circuit are measured with 
and without a perturbing dielectric rod of radius r placed along the axis of symmetry. In the present 
case, the objective was to verify equations (7) and (14) with the experimental results as well as test the 
utility of the available computer codes, ARGUS [8] and SOS [9]. 

A circuit with period L and N periods in length is shorted at planes of mirror symmetry 
(Figure 4). 


Figure 4 - Disk-loaded structure with 5 periods shorted at the ends. 

The resonant frequencies u n are noted at points of maximum signal transmission. The resonant 
frequencies correspond to multiples of half wavelengths along the length of the circuit. If NL = n\ g /2 
and /? = 2x/A ff then 

^ = ( 18 ) 

where n is an integer from 0 to N. For a forward fundamental structure, the highest resonant 
frequency corresponds to the largest 0 and additional periods added to the structure result in lower 
resonant frequencies. 


The interaction impedance is obtained by measuring the dispersion characteristics with a 
perturbing dielectric rod having a small radius introduced along the symmetry axis of the structure. 
The total time average stored field energy in the circuit of length NL is U = (\/2)N LW where W is 
the time average stored energy per period. 

When a perturbing object with polarization P and magnetization M is introduced along the 
axial circuit length, the resulting relative change in resonant frequency / is 


Sf_ 
f ~ 


j (?* ■ E + M* • ll) dv 

V 

2NLW 


(19) 


where the integration is performed over the volume v of the perturbing object. For a thin dielectric rod 
of radius r positioned along the z-axis of the circuit, only the z-component of the electric field should 
be significantly altered. 


^ — ♦ 

In the case of a dielectric rod of relative permittivity e r , the polarization P = (e r — and 

the magnetization M = 0. Therefore the relative change in frequency due to the insertion of the 
dielectric rod becomes 


where 


H 

f 


( ( r ~ VjfQ 

2 NLW 


7T r 


NL 



d *=-5 


-|*»- 2 (€ r -l)c 0 - 


W 


( 20 ) 


NL 

0 

5 


(21) 




The average interaction impedance is then found to be 


7 1 if 1 

0 2 v g f (f r - l)f 0 irr 2 ' 


( 22 ) 


The group velocities were calculated from the derivative of the equation describing the 
dispersion characteristics according to the following relation: 


v g = 2L 


d / 

d 


(23) 


Uniform brass disks with center holes of 0.09375 inch radius and thicknesses of 0.003, 0.005, 
0.008, 0.010, and 0.020 inch were machined. Copper spacers, 0.080 inch thick with center holes of 0.3 
inch radius to act as the cylindrical wave guide, were then stacked with the brass disks to form disk- 
loaded cylindrical wave guide structures. The slot length between disks was 0.16 inch. The lengths of 
the structures measured were from one to ten periods. The resonant frequencies of the structures were 
measured using an HP8510C automatic network analyzer (Figures 5, 6, 7, 8, and 9). The signal 
coupling was achieved by using coaxial probes on both ends of the structure so as to launch a TM 01 
mode. 
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Figure 5 - Transmission measurement using 0.003 inch thick disks. 
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Figure 6 - Transmission measurement using 0.005 inch thick disks. 
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mm (without dielectric tod) ^ (with dielectric rod) 



Figure 7 - Transmission measurement using 0.008 inch thick disks. 
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Figure 8 - Transmission measurement using 0.010 inch thick disks. 


^ iSitiUSl? (without dielectric rod) ^ (with dielectric rod) 



Figure 9 - Transmission measurement using 0.020 inch thick disks. 
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The dispersion characteristics were fitted to a cubic equation in /?£/x (Figure 10). 



Figure 10 - Experimental dispersion characteristics of the disk-loaded structures. 

To determine the interaction impedance, an alumina (e r « 9.2) rod of radius r = 0.0155 inches 
was placed along the length of the symmetry suds of the structure and the resonant frequencies were 
again measured and identified with the corresponding 0 as given by equation (18). 

From the resulting shifts in frequencies and calculated group velocities, the average interaction 
impedances were found from equation (22) and fitted using a cubic spline (Figure 11). 



Figure 11 - Experimental interaction impedances of the disk-loaded structures. 
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COMPUTATIONAL EVALUATION 


Numerical calculations of the resonant frequencies for the 0.010 inch thick disk-loaded wave 
guides were carried out using the two computer codes ARGUS [8] and SOS [9]. These dispersion 
characteristics and average interaction impedances were obtained using the same equations that were 
used for calculation of the experimental results. 

ARGUS is a three dimensional, time dependent or frequency domain electrodynamic and 
plasma simulation code. It is a library of programs built around a unifying framework which utilizes 
the finite difference cell. It contains iterative and non-iterative (direct) Poisson solvers for electrostatic 
problems and also contains an explicit leapfrog Maxwell solver for electromagnetic problems. The code 
can simulate plasmas via an implementation of the particle-in-cell (PIC) method. 

The Self-Optimized-Sector code, SOS, is another three-dimensional, finite difference code for 
simulating electromagnetic processes which involve space charge and electromagnetic fields and their 
interactions. A full set of Maxwell’s time-dependent equations are solved to obtain electromagnetic 
fields and the Lorentz force equation is solved to obtain relative particle trajectories, providing current 
and charge densities. 

A ten period disk-loaded wave guide structure was modeled using the ARGUS computer code 
on the NASA Lewis CRAY Y-MP (Figure 12). 



Figure 12 - Ten disk ARGUS mesh showing the simulated al umina rod through the axis. 

The structure was modeled in a Cartesian coordinate system since the frequency domain 
algorithm in ARGUS can only be used with this geometry. Therefore the cylindrical guide and holes in 
the disks can only be approximated by rectangular elements. The ARGUS code, however, will allow 
the simulated dielectric rod to have a circular cross section with the same radius as that of the actual 
dielectric rod used for the frequency measurements. 

A uniform mesh with a reasonable grid size was chosen to simplify the input file for the code. 
Resonant frequencies were found for the structure with and without a simulated alumina dielectric rod 
(r = 0.0155 inches, e r — 9.2) along the symmetry axis. To obtain more resonances, the number of 
periods in the structure must be increased. Since the lowest resonance frequencies are the frequencies 
found first, lengthening the structure insures that these are the axial resonances. 

The SOS code was run on the NASA Lewis Scientific VAX cluster using one and two periods of 
the structures. For purposes of comparison to ARGUS, a Cartesian geometry was chosen to model the 
disk-loaded cylindrical structure even though the frequency domain algorithm in SOS will work with 
cylindrical geometry. The SOS code used with a Cartesian coordinate system does not allow the 
modeling of the dielectric rod with a circular cross-section, therefore a dielectric rod with a square 
cross-section was used. A uniform mesh was chosen so that the wave guide diameter was an integral 
number of elements such that the area of an element was approximately that of the dielectric rod. 
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A simulated dielectric rod having a rectangular cross section 15% larger than the actual rod 
was used. The difference in the dielectric rod cross-sectional area would be reflected in a difference in 
frequency shift and only have a negligible effect on the values of interaction impedances obtained. By 
using a combination of conductive walls and mirror symmetry boundary conditions at the ends of the 
structures the 0, tt/4, jt/2, 3tt/ 4, and ir points of the dispersion characteristics were found. 


COMPARISON OF RESULTS 

Results obtained through theoretical calculations are compared with those calculated from the 
measured values of resonant frequencies of unperturbed and perturbed 0.010” thick disk-loaded 
structure and those of the two computer codes. The dispersion characteristics are shown for the disk- 
loaded circular wave guide structure obtained using theoretical, experimental, and computational 
techniques (Figure 13). Each set of frequency data was again fitted to a cubic equation in /?X/x. 



0L/n 


Figure 13 - Dispersion characteristics of the 0.010 inch thick disk-loaded structures. 

Resonance frequencies found using SOS are close to those obtained experimentally especially 
near /?X/x = 0.5. The maximum deviation was only 1.63%. Frequency values generated by ARGUS 
were typically lower deviating from 2.3% at pL/it = 0 to 0.75% at /?X/x = 1. Frequency values 
predicted by the theoretical approach agree very well at lower frequencies with measured results but 
tend to be lower at the higher frequencies with as much as a 1.7% deviation. They seem to agree very 
closely to the ARGUS values near /?X/t r = 0.5. 

The average interaction impedances are also shown for the disk-loaded circular wave guide 
structure obtained using theoretical, experimental, and computational techniques (Figure 14). The four 
data sets were fitted to a cubic spline in f}L/n. 

The values of interaction impedance obtained through the use of the SOS code closely agree 
with the theoretical values. This may be due to the similarity in form of the respective dispersion 
characteristics. Likewise, the values of interaction impedance using data from the ARGUS code agree 
extremely well with the experimental results. In this case the form of the respective dispersion 
characteristics are also similar. The similarity in characteristics result in group velocities and hence 
interaction impedances that are nearly the same. 
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Figure 14 - Interaction impedance characteristics of the 0.010 inch thick disk-loaded structures. 


Due to the simplifications in the derivation of equations (14) and (22) the values of interaction 
impedance in all cases are very high and serve only as an index of comparison between the different 
approaches. Also the method of calculating the group velocities lends uncertainty to the values of 
interaction impedances. In general, the behavior of impedance is consistently similar in all the four 
cases. The minor deviations of the computer values from the experimental values can be attributed to 
the computational approaches used in ARGUS and SOS. 


The theoretical results incorporated several assumptions in the theory. The Helds in the slot 
are assumed to be radial TEM and for X g t -f- s this assumption is justified. Space harmonics which 
would introduce a sin(/? n s)/(/? n «) factor for the nth harmonic into the analysis have been ignored. 
Admittance matching across the slot at p = a is done in the average sense using the space factor tj to 
account for the electric field being zero on the disk edge. Actually a point by point match cannot be 
done and so inaccuracies will inevitably occur in these results. 


The design of structures using the theoretical approach in association with an appropriate 
computer code can be greatly facilitated and can provide an insight in the performance of a structure 
prior to its being fabricated and tested. The disk-loaded structure is a good candidate for use in 
miniature traveling-wave tubes. With the proper choice of dimensions a disk-loaded device with good 
interaction impedance operating at a nominal beam voltage may be possible. 


CIRCUIT MODELING 

The disk-loaded structure lends itself readily to a circuit modeling approach. A way to 
describe one period of the structure is to use a series of four cylindrical wave guides, two cylinders with 
the radii equal to the radius a of the disk hole and having a length equal to the half the thickness f/2 
of the disk between two cylinders each having radii equal to the radius b of the disk-loaded structure 
and having lengths equal to half the gap spacing s/2. 

Harrington [10] provides equivalent circuits of a cylinder for both TE and TM modes. In slow- 
wave structures for traveling-wave tube applications we are interested in the TM modes only since 
these provide a non-zero axial electric field which can couple to an electron beam. 
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For the lowest order cylindrical TM mode, a TM 01 mode, the circuit equivalent values per unit 
length for impedance Z and admittance Y for a cylinder of radius r are given by 

k 2 

Z = jX = ju>n - j -fe (24) 

and 

Y = jB — jwe (25) 


where k c = 2.405/r is the cutoff wave number. 

The values for the characteristic propagation constant 7 and characteristic impedance Z 0 for 
the TM 01 equivalent circuit can also be found using general transmission theory. 

The characteristic propagation constant and impedance are given as 

7r M = <ZY = jP = jH 1 ~(k c /k) 2 (26) 

and 

ZoTM = WlY = Ul = 4^1 - (*c/*) 2 (27) 


where k = w-^juc is the wave number. 

Circuit elements having values corresponding to the 0.010 inch thick disk-loaded wave guide 
structure were used to model cylinder sections. The circuit modeling was performed using Touchstone. 
Touchstone is a linear, frequency domain circuit program which can perform simulation and 
optimization of an electrical circuit. The program resides on a SUN SPARC workstation at NASA 
Lewis Research Center. 

Circuit elements were also needed to account for discontinuities at the ends and between 
dissimilar cylinders. Circuit elements with a variable series LC impedance and a variable parallel LC 
admittance were chosen to model the discontinuities. The entire circuit was kept symmetric with 
respect to the center of the structure (Figure 15). 
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Figure 15 - Touchstone model of the 0.010 inch disk-loaded circular wave guide. 

The values for the inductances and capacitances for the variable elements were found by 
optimizing the circuit such that maximum transmission occurred at the 0L = 0 and /?£ = x 
experimentally measured cutoff frequencies of the structure (Figure 16). The transmission is lossless at 
these frequencies since no resistors were introduced into the circuit. 


12 




EEsof - Touchstone 
o MAG(Sa] 



Figure 16 - Tr ansmis sion simulation using Touchstone. 

This effort shows that it is possible to model a period of a given slow-wave structure using the 
0L = 0 and (3L = tt known resonance frequencies. It is a matter of adding resistances to account for 
the losses observed. 
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APPENDIX A 


Math cad program to calculate frequency and interaction impedance for a disk-loaded circular wave guide 
TM 01 mode given the inner radius, outer radius, slot spacing, disk thickness, and axial wavenumber. 


inner radius 

a := .09375- .0254 

speed of light 

c := 2.99792- 10 8 

TE 11 cutoff 

outer radius 

b := .3-. 0254 

electronic charge 

q ;= 1.60218-10' 19 

1841 in - 9 
c -10 

slot spacing 

s := 2- .08- .0254 

electronic mass 

m:r 9.1095 10 ' 31 

2 - *b 

TM 01 cutoff 

disk thickness 

t := .020- .0254 

permeability 

p := 4-X-10 ' 7 

2405 in 9 
c -10 

period 

d := S + t 

permittivity 

s := 8.854215- 10' 12 

2 -x-b 





TE 21 cutoff 


t 


j r 

3-054 9 

spatial factor 

q := - 

axial wavenumber 

P = (. 8 )-- 

c -10 

d 


d 

2 -x-b 


Define Bessel function combinations 

Z00(k,a) = JO(k-a)- YO(k-b) - Y0(k a) JO(kb) ZOl(k.a) ;= JO(k-a)- Yl(k-b) - YO(k-a) Jl(k-b) 

Z10(k,a) = Jl(k-a)-Y0(k-b) - Y1 (k a) JO(k b) Zll(k.a) = Jl(k a) Yl(k-b) - Y1 (k-a)-Jl(k b) 

Define Bessel function combination derivatives 


ZOO'(k.a) ;= -ZlO(k.a) - -Z01(k,a) 

a 

Define a and the derivative a‘ 

a(k>a):= J_.zio(M) 

k-a Z00(k,a) 

Define modified Bessel functions 


ZlO'(k.a) := Z00(k,a) - --ZlOCk.a) - - Zll(k,a) 
k-a a 


, „ N / ZlO'(k.a) ZOO'(k.a) 1 

a'(k,a) := a(k,a)- - ■ 

\Z10(k,a) Z00(k,a) k-a 


j = 0.. 15 


I 0 (x) := ^ 


2 j 


2-j+l 


J 


(J) 


i 1( x):,yA 2 i 


I2(x) := I0(x) - ~Il(x) 
x 


Define modified Bessel function derivatives 

I0'(x) := II (x) 

Define 4> and the derivative <t>' 

■ ^ . 1 Il(x-a) 

4>(T,a) := 

x-a IO(x-a) 

Solve for k using $=a 

k := 300 GIVEN 


Il'(x) := I0(x) - 


II (x) 


<t>'(T.a) := 4>(t»a)- 




4>WP - k ,a/- 


. a(k,a) 
i - n 


Il’(T-a) IO'(Ta) 

II (x a) IO(x-a) x-a 


k := FIND(k) 


Calculate frequency, x, a, and 4> 


f := — — 
2 -x 

x:= Jp 2 -k 2 

<t>(x,a) =0.441851 

f- 10 ' 9 =15.693878 

x =440.44614 

a(k,a) =0.392757 
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Calculate the phase velocity 

V : = c~ - =0.598352 

P c 


Calculate the beam voltage 


V = 


1 


2 

1 - 1 - 

2 

c 


m 2 

— c 


= 126767.243147 


Calculate the group velocity 


Calculate the power 




vg := 


Pa P-a 1-n 


— =0.023221 


P :■ iti. 1 - 

2 XT 


— (Il(x-a) Il(x a) - I0(xa)-I2(ta)) 
P 


dpdk c 

Calculate the interaction impedance 


P =3.635789-10 


dielectric rod radius r = 0.0155-.0254 Ez2 = I0(xr)I0(x-r) - Il(t r) Il(x r) Z = 


Ez2 

2ppP 


Z =458.531069 


Calculate the matching impedance 

zi ■, c .i 2 oo. 10 < T :;i^i 

k x-all(xa) 


— 9 10 11 =19.188538 
Zi 


Calculate the lower cutoff frequency due to the outer radius b 

k ;= 300 GIVEN J0(kb)=0 klow = Find(k) flow = klow-10' 9 — flow = 15.058052 

2jc 

Calculate the upper cutoff frequency due to the inner radius a 

GIVEN Jl(k a) Y0(k b) - J0(k b) Yl(k a)-0 khigh = Find(k) fliigha = khigh 10‘ 9 — fhigha = 17.65171 

2* 

Calculate the upper cutoff frequency due to the period d 

fhighd := - 10‘ 9 fhighd = 32.785652 


Solve for the dispersion characteristics 

n = 10 l=0..n P(l) = -■ - k(l) = root 

n d 


a(Re(k),a)-(l-n) 


M 


p(l) 2 -Re(k) 2 ,a, 



Plot of the dispersion characteristics 


1 f(l) 10 9 


0 


15.306505 

1 


15.314177 

2 


15.336804 

3 


15.373274 

4 


15.421881 

5 


15.480517 

6 


15.546883 

7 


15.618684 

8 


15.693789 

9 


15.770333 

10 


15.846771 



f(l) 


ck(l) 

2 % 
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APPENDIX B 


Mat head program to calculate frequency and interaction impedance for a disk-loaded circular wave guide 
TM fll mode given the inner radius, outer radius, slot spacing, disk thickness, voltage, and current. 


inner radius 

a := .09375-. 0254 

speed of light 

c 2.99792- 10® 

TE 11 cutoff 

outer radius 

b = .3-. 0254 

electronic charge 

q := 1.60218 10* 19 

1.841 , 9 ,, 

c-10 =11.527602 





2 -x-b 

slot spacing 

s := 2- .08- .0254 

electronic mass 

m := 9.1095-10 31 

TM 01 cutoff 

disk thickness 

t := .010- .0254 

permeability 

p := 4-x-lO ' 7 

2.405 9 

c- 10 =15.059144 

period 

d := S + t 

permittivity 

e := 8.854215- 10' 12 

2-*-b 





TE 21 cutoff 

spatial factor 

t 

” = ~ d 

beam voltage 

V := 60000 

3.054 9 

-c-lO'* - 19.122921 



beam current 

I:= .001 

2 -x-b 


Calculate the phase velocity 







llsl) 2 



2\ 

2 

\ m 'C / 


1 

me 

- 1 

2 

.J 

lq-vi 


Define Bessel function combinations 


- =0.457026 
c 


ZOO(k.a) := JO(k-a)- YO(k-b) - Y0(k-a)-J0(k b) Z01(k,a) := JO(k-a) Yl(k-b) - Y0(k-a)-Jl(k-b) 

Z10(k,a) := Jl(k-a)- YO(k-b) - Yl(ka)-J0(k-b) ZU(k.a) := Jl(k-a)- Yl(k-b) - Yl(k-a)-Jl(k-b) 


Define Bessel function combination derivatives 


Z00'(k,a) := -ZlO(k.a) - -Z01(k,a) 
a 

Define a and the derivative a' 

a(k,a) = _* — ZlOClc^a) 
ka Z00(k,a) 

Define modified Bessel functions 


j = 0. 15 


I0(x)=5] 


2j 


J 


r or 


Define modified Bessel function derivatives 

IO’(x) := II (X) 

Define and the derivative <t>’ 

1 Il(x-a) 


4>(T,a) := 


x a IO(i-a) 


ZlO’(k.a) := ZOO(k.a) - - -ZlOCk.a) - -Zll(k,a) 
k a a 


a'(k,a) := a(k,a)- 


I ZlO'(k.a) ZOO’(k.a) 1 


ii ( x),VA?L. 

. j !'0 + 


\Z10(k,a) Z00(k,a) ka 

2 j + 1 


I2(x) = I0(x) - - II (x) 

X 


1 )! 


Il’(x) := I0(x) - 


Il(x) 


f(x,a) := <t>(x,a)- 


Il'(T-a) IO'(xa) 
Il(t-a) I0(ta) xa 


Solve for k using 4>=oc 
k = 300 GIVEN 



a(k.a) 

1 - T] 


k := FIND(k) 
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p=-k 

V 


Calculate frequency, p, x, a, and $ 

f ■= — — 

2 * 

f- 10 9 = 15.79377 p =724.277576 

Calculate the group velocity 

dpdk= ^ + Taa '< k - a > 1 


X := Jl P 2 - 1 


x =644.211229 

Calculate the power 


a(k,a) =0.37084 
4>(x,a) =0.394017 


vg := 


P a p-a f (x,a) 1 - n 


— =0.023943 


P , 

2 x-x 


— (Il(x a)-Il(x a) - I0(x a) I2(x a)) 
P 


P =5.91736*10 


dpdk c 

Calculate the interaction impedance 

dielectric rod radius r := .0155- .0254 Ez2 := IO(x-r) IO(x-r) - Il(x r) Il(x r) 

Calculate the matching impedance 


Z = 


Ez2 

2-p-P-P 


Z = 163.687697 


P ^ I0(x a) - 1 

Zi := c - 200 — - — 

k 


— -9- 10 11 =15.705308 
Zi 


xa-Il(xa) 

Calculate the lower cutoff frequency due to the outer radius b 

k = 300 GIVEN J0(kb)=0 klow = Find(k) 

Calculate the upper cutoff frequency due to the inner radius a 


flow := klow- 10' 9 — 
2 k 


-9 c 


flow = 15.058052 


GIVEN Jl(k-a)-Y0(k-b) - J0(k b) Yl(k a)-0 khigh = Find(k) fhigha = khigh 10 — fhigha = 17.65171 

2-s 

Calculate the upper cutoff frequency due to the period d 

fhighd = 34.71422 


fhighd ;= — — • 10' 9 
2-d 


Solve for the dispersion characteristics 

1 K 


1 = 0.. n P(l) := k(l) = root 

n d 


n := 10 

Plot of the dispersion characteristics 

400 


«(R e (k),a) - (1 - r, ) ■ 4» W P ( 1) _ Re(k) ,a/,kj f(l) = 


c-k(l) 
2 k 


f(l) 10 


-9 




0 


15.185202 


1 


15.19403 


2 


15.220011 


3 


15.261706 


4 


15.321726 


5 


15.38878 


6 


15.4641 14 


7 


15.544958 


8 


15.628804 


9 


15.713518 

X 

10 


15.797392 
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APPENDIX C 


Mathcad program to calculate the inner radius and interaction impedance for a disk-loaded circular 
wave guide TM 01 mode given the frequency, outer radius, period, disk thickness, voltage, and current. 


outer radius 

b := .006625 

speed of light 

c := 2.99792- 10® 

TE 11 cutoff 

beam voltage 

V := 5000 

electronic charge 

q := 1.60218 10' 19 
.31 

1 84 1 o 

-c-10" 9 -13.258918 

2-x-b 

beam current 

I .001 

electronic mass 

m:= 9.1095 10 

TM 01 cutoff 

frequency 

f:* 30- 10 9 

permeability 

p = 4-Jt 10‘ 7 

^^-c-10' 9 = 17.320857 

disk thickness 

t := .00006858 
2.405 

permittivity 

s= 8.854215 10 12 
2-x-f 

2-x-b 

TE 21 cutoff 

3.054 .9 

minimum outer radius c 

2-x-f 

= 0.003825 

k:= 

c 

c-10 * -21.994967 

2-x-b 


Calculate the phase velocity, k, p, and x 


/ _ ..\2 




2_q-V 


/ n\ 

2 

\ m-c / 


( m-c 2 ' 

1 + 

\q-vi 

- 1 

2 

[J 


- =0.139208 
c 


p := - k X J p 2 - k 2 

v 


Choose the period d 

maximum period — =0.000696 period d = .4 — d =0.000278 

P P 

Define Bessel function combinations 


spatial factor 



tj =0.246493 


ZOO(k.a) := J0(k a)-YO(kb) - Y0(k-a)-J0(k-b) Z01(k,a) := JO(k-a)- Yl(k-b) - Y0(k-a)-Jl(k-b) 

Z10(k,a) := Jl(k-a)-Y0(k-b) - Y1 (k-a)-JO(k-b) Zll(k,a) := Jl(k-a)- Yl(k-b) - Yl(k a) Jl(k-b) 


Define Bessel function combination derivatives 


ZOO'(k.a) := -Z10(k,a) - --Z01 (k.a) 

a 


ZlO'(k.a) := ZOO(k.a) - --Z10(k,a) - - Zll(k.a) 
k a a 


Define a and the derivative a' 


a(k,a) 


_1_ Z10(k,a) 
k a Z00(k,a) 


Define modified Bessel functions 


a'(k,a) .= a(k,a) 


/ Z10'(k,a) 

\Z10(k,a) 


ZOO'Ck.a) 

Z00(k,a) 



j = 0.. 15 I0(x):=^] — 


2j 


J 


j!j!-2 


2j 


noo-S— 


2j + 1 


j! (j + 1)1 2 


2‘j + 1 


I2(x) := I0(x) - -Il(x) 
x 


Define modified Bessel function derivatives 


I0'(X) := II (x) 


Define 4> and the derivative 4>* 

1 Il(x-a) 

<Kx,a) := 

x-a IO(x-a) 


Il'(x) := I0(x) - 

X 


f(x,a) = 4»(x,a)- 


/ il'(x-a) 
\Il(x a) 


IO'(xa) 1 \ 
IO(x-a) xa/ 
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Solve for the inner radius a 

a = .5 b GIVEN <t>(x,a)- a(k - - a = FIND(a) 

i - n 

Calculate the group velocity 


Calculate 41 and a 

a =0.004015 <|»(T,a) =0.054157 

Calculate the power 


a(k,a) =0.040808 


dpdk = 


k-a x-a a'(k,a) 1 

p a P-a 4>'(t,a) 1 - r\ 


P = 


------a 2 - - (Il(x a) ll(x a) - I 0 (x a) I 2 (x a)) 

2 T x 


vg = 


c 

dPdk 


^ =0.007476 
c 


P = 18288.69785 


Calculate the interaction impedance 

dielectric rod radius r- .01 55- .0254 Ez2 = I0(x r) I0(x r) - Il(x r) Il(x r) 

Calculate the matching impedance 


Ez2 

Z := 

2-p-p-P 


Z =2.881309-10 


Zi.o^OO 10 ^- 1 
k xall(xa) 

Calculate the Pierce parameter 

.33333 


C = [-- 

4 V 


C =5.243035-10 


-7 


- 9- 10 11 =36.490534 
Zi 

Calculate the rf power output 

2CV-I =0.000005 


Calculate the lower cutoff frequency due to the outer radius b 

k = 300 GIVEN J0(kb)=0 klow = Find(k) 

Calculate the upper cutoff frequency due to the inner radius a 

GIVEN Jl(k a) Y0(k b) - J0(k-b)- Yl(k a)-0 khigh = Find(k) 

Calculate the upper cutoff frequency due to the period d 


flow = ldow- 10' 9 -— 
2 * 


fhigha = khigh- 10’ 9 — 

2-Jt 


flow = 17.3196 


fhigha =31.665994 


fhighd = - IQ 9 fhighd = 538.762802 
2-d 

Solve for the disperion characteristics 

n := 10 1 = 0..n 

Plot of the dispersion characteristics 


P(l) = k(l) = root 
n d 


a(Re(k),a) - (1 - n) 4>WP(1) ' Re(k) 


a 



f(l) 10 


-9 



0 


18.646066 


1 


25.926408 


2 


28.514995 

_ 

3 


29.487839 


4 


29.998627 


5 


30.300609 


6 


30.469927 


7 


30.556855 


8 


30.601727 


9 


30.62655 

t 

10 


30.641434 


l ,a) ,k 


f(l) 


c-k(l) 

2 % 
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APPENDIX D 


Mathcad program to calculate the interaction impedance for a disk-loaded circular wave guide 
TM 01 mode given experimental or computational data. 


inner radius 

a := .09375-. 0254 

speed of light 

c := 2.99792- 10* 

Read frequency data file 

outer radius 

b = .3-. 0254 

electronic charge 

q:= 1.60218 10 19 

M := READPRN ( dcwgOO 10) 

disk thickness 

t := .003. 0254 

electronic mass 

m := 9.1095 10 31 

_ _.7 

Unperturbed frequencies 

f^M^-10 9 

slot spacing 

s = .16- .0254 

permeability 

^ := 4 x10 

* 

period 

1 := S -t- 1 

permittivity 

e:= 8.854215- 10 12 

f=f 

Perturbed frequencies 

dielectric rod radius 

r := .0155-. 0254 

dielectric constant 

er := 9.2 

fd=M <1> 10 9 

Frequency shifts 8f(i) := f f - fd 

Least squares cubic fit in p 

number of points 

N = length(f) 

fd := fd 

i=0..N-l D:=3 d 

Calculate the coefficients 

:=0..D P r 

= i+ p=p 

N 1 



to := 2-X-f 


coeff := (x T x) 1 -(x T -m) Q(x) := ^coeff d x d Q'(x) := J]ooeCf d -d-x d " 1 


Calculate the phase velocity, group velocity, and equivalent beam voltage 

«(Pi) 

vp(i) : = V ' 


Pi 


Vg(i) := a Pj 


V(v) := 


V 1 


1- 


c 2 / 


- 1 


2 

c 

m — 

q 


Calculate the interaction impedance and gap factor Plot of frequency, fit, and perturbed frequency 


8f(i) 


Z(i) := 


2 2 
•vg(i) (er- l) e x r 


M(i) := 



l A 

sin 

Pj— 


[ ' 2/ 




Values of frequency, fit, and interaction impedance 


i + 

1 

f; 10 9 


Z(i) 

i 


15.178125 


15.175524 


58614.570394 

2 


15.228125 


15.232496 


6627.354587 

3 


15.325 


15.321675 


1942.569347 

4 


15.425 


15.432365 


726.471087 

5 


15.55625 


15.553868 


363.295224 

6 


15.684375 


15.675488 


216.833099 

7 


15.784375 


15.786528 


124.327911 

8 


15.875 


15.876289 


111.09255 

9 


15.928125 


15.934076 


139.059421 

10 


15.953125 


15.949191 


-344.364848 
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